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Abstract. This review describes the evidence for small-scale structure, 'clumping', 
in the radiation line-driven winds of hot, massive stars. In particular, we focus on ex- 
amining to what extent simulations of the strong instability inherent to line-driving can 
explain the multitude of observational evidence for wind clumping, as well as on how 
to properly account for extensive structures in density and velocity when interpreting 
the various wind diagnostics used to derive mass-loss rates. 



1. Introduction 

Hot, massive stars possess strong stellar winds driven by line scattering of the star's 
intense continuum radiation field. The first quantitative description of such line driving 
was given in the seminal paper bv lCastoretaD (119751 . 'CAK'), assummg a stationary, 



homogeneous, and spherically symmetric wind. But despite the co nsiderable s uccess 
of subsequent work assuming such a smooth, steady outflow (e.g.. IVink et al.ll200Q) . 



the theoretical as well as observational evidence for a time-dependent, inhomogeneous 
wind is irrefutable (iPuls et al.ll2008l : IHamann et al.ll2008l) . This review summarizes the 
current status concerning small-scale inhomogeneities, 'clumping', in the winds of sin- 
gle hot stars; for an overview of large-scale wind inhomogeneities, induced by for ex- 
ample rapid rotation or a strong magnetic field, see lPuls et al.l (I2008f) . After first outlin- 
ing the theoretical basis in terms of a strong instability inherent to line driving, we next 
concentrate on confronting predictions of corresponding instability models with a sub- 
set of the large body of direct and indirect observational evidence for wind clumping. 
In particular we focus on the consequences of extensive density and velocity structure 
when interpreting various mass-loss diagnostics. 



2. Theory: the line driven instability 



Linear stability analysis showed already early on that the line-driven winds of hot stars 
should be unstable for velocity perturbatio ns on scales near and below the Sobolev 
length Lsob = Vib/idv/dr) (MacGregor e t aI]ll979l : ICarlberglll980i : IOwocki & Rvbicki 
[l984) . The operation of this strong, intrinsic 'line-driven instability' (L DI) has since 
been c onfirmed by direct numerical modeling of the time-dependent wind dOwocki et al. 
19881 : iFeldmeieil 19951 : iDessart & Qwockil 2005al) . Such simulations typically show 



that the non-linear growth of the LDI leads to high-speed rarefactions that steepen into 
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strong reverse shocks, whereby most of the wind material is compressed into dense and 
spatially narrow 'clumps' that are separated by large regions of much lower densities. 
This characteristic structure is the theoretical basis for our current understanding and 
interpretation of wind clumping. 

The left panel of Fig. [T] illustrates this typical structure by plotting density and 
velocity snapsho ts of a 1-D LDI simul ation computed from an init ial steady CAK 
model following iFeldmeier et all (1997) (see ISundqvist et all 1201 Ibl for model spe- 
cific details). T he line force is calculated using the 'Smooth Source function' method 
(SSF. Owocki & Puis 1996) . which allows one to follow the non-linear evolution of the 
strong, intrinsic instability, while als o accounting for the stabilizing Une drag effect of 
the scattered, diffuse radiation field dLucyl [19841) . The simulation h ere further intro- 
duces base perturbations from Langevin turbulence (see IFeldmeier et al.i,1997.) . which 
induce structure somewhat close r to the wind base as cornpared to simulations with 
self-excited structure formation ( Runacres & Owockill2002l. 12005 ). But despite these 
seed perturbations, the line drag from the SSF force greatly reduces, even eliminates, 
the instability near the base, thus allowing the lowermost wind to remain smooth and 
steady (Fig.[T||. 




r/R, n/v. 



Figure 1. Left: Snapshot of density and velocity structure in a 1-D LDI simula- 
tion (blue), as compared to the smooth CAK start model (red, dashed), see text. Mid- 
dle: Observed lUE spectra of the C iv /1/11548, 1551 resonance doublet in ^ Pup, as 
function of line-of-sight velocity of the blue component. Overlayed are profiles cal- 
culated from CAK and LDI models, where the latter synthetic spectr u m is compute d 
using an extension of the technique developed in ISundqvist et al.l (1201 Ot 1201 Ibh . 
which takes 1-D snapshots and phase them randomly in patches of a parametrized 
angular size, here of 1°. Right: Density contours from a 2-D LDI simulation, in 
which the A<p -12° wedges repr esent a clockwise time seq uence starting from the 
initial CAK model. Adapted from lDessart & OwockH (120031) . 



The presence of strong, embedded wind shocks in LDI simulations provides a 
generally accepted explanation for the X-rays observed from sin gle hot, massive stars 
without very strong magnetic fields, as further discussed by, e.g. lOwocki et all (|201 1 ) 
in this volume. In simulations with self-excited structure though, there is not enough 
material going through strong enough shocks to produce the observed amount of X- 
rays. But as shown by Feldmeier et all (|l997h . introducing base perturbations such as 
those considered here enhances the velocity dispersion of the wind clumps, whereby 
clump-clump collisions create regions of relatively dense, hot gas. This then increases 
the X-ray emission to levels more comparable to those typically observed. 

Because of the computational expense of calculating the line-force at each time 
step, LDI simulations have generally been limited to 1-D. More realistically though, 
thin shell instabilities and associated effects can be expected to break up the spherical 
shell structure into a more complex multi-dimensional structure. First attempts to con- 
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struct 2-D LDI simulations have been carried out bv lDessart & Qwockil (I2003L 12005 ah . 
Such models typically show that the LDI first manifests itself as strong density com- 
pressions mimicking corresponding 1-D simulations. But as these initial shell structures 
are accelerated outwards, they become disrupted by Rayleigh-Tay lor and thin-shell in- 
stabilities, which in the simulations by Dess art & Qwockil (l2003h operate all the way 
down to the azimuthal grid scale (Fig. |2j right panel). Moreover, these 2-D simula- 
tions have higher velocity dispersion than corresponding 1-D ones. Similarly as for the 
models with base perturbations just discussed then, this suggests an increase also in the 
amount of X-rays produced. 

However, it may well be that these initial 2-D models exaggerate somewhat the 
level of lateral disruption, as they do not yet properly treat the lateral component of the 
diffuse line drag. Presuming this could damp azimuthal velocity p erturbations on scales 
below the lateral Sobolev length = rvth/vr (iRybicki et al.lll990h . lateral breakup may 
be prohibited below scales of order A(p ~ L^/r ~ Vth/vp ~ 0.5°. Future work should 
examine this issue by an adequate incorporation of the lateral line-force into multi-D 
LDI simulations. 



3. Observations 



3.1. Line profile variability 



The most direct evidence for clumping in hot star winds comes from observations of 
narrow, spectra l subpeaks superimposed on b road optical emission lines in Wolf-Rayet 
(WR) stars (Moffa t et al.lll988l : lRober3ll994 . Careful time-monitoring have revealed 



how these subpeaks, which trace local wind density enhancements, systematically mi- 
grate from line center towar d line edges on tirne scales that can be associated with the 
wind acceleration. Indeed, iDessart & Owockj (l2005b ) showed that synthetic spectra 



calculated directly from LDI models can match well the observed migration in WR 
winds; however, the inferred acceleration scale typically predicts either a more ex- 
tended acceleration zone than dynamically expected for line driven winds, or a much 
larger core radii than expected for WR stars (see also'Lepine & Mo fFad 19991) . 



While predominate ly observed in dense WR winds, Eversber g et al.l (Il998 ) and 



Lepine & MoffatI (l2008h have provided first evidence of emission substruct ures also in 



0-stars, suggesting the phenomenon is universal in hot star winds. Moreover. Eversbe rg et al. 
(0*998) pointed out that the density enhancements in ^ Pup always seemed to appear 
very near the stellar surface, thus providing an early indication of clumping close to the 
wind base (see Sect. [5]l. Further monitoring of such emission lines formed mainly in 
the inner wind would provide crucial information regarding clumping and wind accel- 
eration in this key region. 

3.2. The absorption troughs of saturated P Cygni lines 

Instability simulations predict extensive structure not only in density, but also in veloc- 
ity (Fig.[T] left panel). The most clear-cut evidence for such velocity structure in hot star 
winds comes from the extended regions of zero residual flux observed in saturated UV 
resonance lines (Fig.[T] middle panel). Already Lucy (1982) suggested that the physical 
origin of these 'black troughs' is the non-monotonic wind velocity field, which leads 
to enhanced back-scattering in multiple resonance zones and thereby to a systematic 
reduction of the blue-shifted emission. The middle panel of Fig. \T\ demonstrates that 
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the velocity structure predicted by LDI simulations well reproduces the observed char- 
acteristics, whereas the monotonic velocity law of the CAK comparison model fails 
to reproduce in particular the blue edge of the line profile. This clearly indicates the 
presence of material with velocities well above the wind terminal speed Voo that sets the 
maximum blue-shifted absorption in a CAK modelQ 



4. Diagnostics: deriving mass-loss rates from clumped hot star winds 

Due to the rather steady wind base, the average mass-flux in the simulations described 
in Sect. |3] is actually not much affected by the LDI, which may indicate only second 
order effects on theoretical mass-loss rates predicted by steady-state models. But let 
us now discuss, arguably, the key question concerning wind clumping, namely how the 
structures in density and velocity affect the various diagnostics used to infer mass-loss 
rates from observations. 



4.1. Optically thin clumping 

Traditionally most diagnostic studies aiming to derive mass-loss rates have assumed 
that the wind consists of statistically distributed optically thin clumps embedded in 
a void background medium, and neglected any disturbances on the velocity field. The 
main effect of such optically thin clumping is that diagnostics having opacities that scale 
with the local density squared, such as H„ and thermal IR/radio free-free emission, are 
stronger than in a smooth wind with the same mass-loss rate. The scaUng invariant is 
yffd^, i-S- mass-loss rates derived from Hq, and smooth wind models are overestimated 
with the square-root of the clumping factor /d = (p^)/(p)^ in the Hq, line forming 
region. In contrast, diagnostics having opacities that scale linearly with density, such 
as UV resonance lines, electron scattering wings, and bound-free absorption of X-rays, 
are not directly affected by optically thin clumping. 

For WR-stars, optically thin clumping has been accounted for since the pioneer- 
ing work by iHillieii (119911) . Simultaneous fitting of electron scattering wings (scaling 
linearly with density) and emission peaks (density-squared) of optical recombination 
lines indic ate fci ~ 4-20, c orresponding to mass-loss reductions by factors of ~ 2-4 (see 
review by ICrowthetll2007h . For O stars, electron scattering wings are too we ak for this 
technique to be of use. But combined H^/IR/radio studies ( Puis et all 12006*) yield up- 



per limit mass-loss rates that, similarly, are factors of ~ 2-4 lower than earlier Hq. rates 
based on smooth wind models, and r ecent X-ray line at t enuati on studies confirm such 



reductions also on an absolute scale (ICohen et al.ll2010l. 1201 1 ). H owever, un saturated 



UV r esonance lines often suggest much more drastic reductions (iBouret et al. 20031, 
I2OO5I : IPrinja et all l2005b . The most prominent exam ple in this respect has been the 
analysis of phosphorus v (Pv) in 40 Galactic O stars by iFuUerton et all (l2006t) . who de- 
rived values of the mean ionization fraction times the mass-loss rate, (^pv)M, that were 
factors of 10 . . . 100 lower than corresponding M values derived from smooth models 



'The only way these features can be modeled within the context of a smooth outflow is by introducing a 
highly supersonic 'micro-turbulence' ( lHamannlll98lh . typically on the order of ~5-10% of the terminal 
speed. 

^They can be indirectly affected though, th rough a modified wind ionization balance due to increased 
recombination rates, e.g. iBouret et all d2003l) . 
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and H„/radio emission. Since state-of-the-art non-LTE model atmospheres that include 
optically thin clumping typically predict <^pv) « 1 for early to mid type O stars, these 
Pv results seem to imply extremely low mass-loss rates, thus challenging the validity 
of the line-driven wind theory. 

4.2. Optically thick clumping: porosity and vorosity 

It is important to realize that the above results all rely on the assertion that individual 
clumps are optically thin; if this is not true for the investigated process, additional ef- 
fects become important in the radiation transfer. Namely, optically thick clumps lead to 
a local self-shielding of opacity within the clumps, which in turn allows for increased 
escape of radiation through porous channels in between the clumps. The essential effect 
of such porosity is to reduce the medium's effective opacity as compared to an optically 
thin clump model (Feld meier e t al. 2003 ; Owocki et al. 200d); Fig. |2] graphically illus- 
trates this general increase in transparency for porous media. The trade-off between 
porosity and mass-loss rate thus opposes the trade-off for density-squared diagnostics 
between optically thin clumping and mass-loss rate , i.e. neglecting porosity in cases 
where it is important can lead to underestimated mass-loss rates. 




Figure 2. Left, middle: Randomly generated isotropic clumps of length scales 
{c\ ~ O.lr and velocity-stretch porosity with hoo - 0.5/?* (left) a nd hca = 5.0/?,. 
middle ), that are lit from behind by a uniform source. Adapted from lSundqvist et al] 
201 lal) . Right: Phenomenological velocity stair-case model vs. wind mass, where 
the dense clumps have velocity spans Sv that are separated by an amount Av. Here 
the associated velocity filling factor is /vei - 6v/ Av - 1/10 and the straight line 
represents a smooth CAK velocity law. Adapted from lOwockil (l2008l) . 



Porosity in continuum absorption of X-rays? For continuum processes, the key pa- 
rameter controlling to which extent porosity might be important is the so-called poros- 
ity length, defined as the mean free path between clumps and also equal to the clump 
length scale times the clumping factor, h - lafd- Simplified porosi ty wind models 
have been applied to bound-free X-ray attenuation in 0-star winds (iFeldmeier et al.l 
20031 : lOskinova et al.ll2004. ,20061 : lOwocki & Cohenll2006l) . The specific results of these 



studies depend on whether isotropic or anisotropic clumps are assumed, where the lat- 
ter tends to enhance porosity effects (lOskinova et al.ll2006l : ISundqvist et aLll2011al) . A 
key result for isotropic porosity is that it is quite difficult to make individual clumps 
optically thick, simply bec ause the total integrated optical depths are quite low for X- 
ray absorption in O stars (ICohen et al.l 12010 , 1201 lb . For 'velocity stretch' models in 
which the porosity length scales with the local wind velocity, h - /Joo(v/voo), terminal 
poro sity lengths well above a stellar radius typ ically are required for a significant ef- 
fect (ICohen et al.ll2008l : ISundqvist et al.ll2011al) . Therefore it seems rather unlikely that 
porosity from LDI structure could have a significant effect on X-ray line profiles; even 
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though velocity stretching and clump-clump collisions in 1-D LDI models can cause 
clump separations in the outer wind to become on the order of a stellar radius, clump 
separations at lower radii tend to be much smaller, on the order of the Sobolev length 
Lsob = Vth/{dv/dr) » (vth/v)/?* « 0.01/?*. The left panel of Fig. |3] demonstrates ex- 
pUcitly that X-ray line profiles with an opacity computed directly from the LDI density 
structures described in Sect. |2] are very similar to profiles computed using a smooth 
CAK model, indicating small porosity effects. 

Vorosity in resonance line absorption? Contrasted to the case of continuum absorp- 
tion discussed thus far, clumps can easily become optically thick in the in herently very 
strong UV resonance lines, as first pointed out by lOskinova et al.l ( 2007 ). In a rapidly 



accelerating stellar wind, each line photon can only interact with the wind material 
within a very narrow spatial range, due to the narrow Doppler width of the line profile. 
This makes it possible for line photons to leak through the wind via 'porous' channels 
in velocity space, without ever interacting with the optically thick clumps (Fig. [21 right 



most panel); hence this effect has been dubbed velocity porosity, or 'vorosity' (lOwocki 



120081) . For a simple model with zero thermal speed, such vorosity may be characterized 
using solely a velocity filling factor /vei = 6v/Av, as defined in Fig. |2] The basic point 
though, is that this photon leakage in velocity space produces the same general effect 
as spatial porosity, namely a reduction in effective opacity as compared to an optically 

thin clump model. 

iPrinja & Massal (1201 Ol) found empirical evidence for such optically thick clumping 
by analyzing the Si iv resonance doublet in B supergiants. As these resonance lines are 
unaffected by optically thin clumps, the optical depth ratio between the blue and red 
line components should reflect the underlying atomic physics, i.e. t^/t^ = fb/fr = 2.0 
where / is the os cillator strength, if clum ps were indeed optically thin. However, the 
ratios inferred by 'Prinja & Massal ( 20101) spread over the range ~ 1 to 2, with a mean 



Tb/Ty = 1.5 ± 0.3 derived from their core sample of 25 stars. This result is a clear 
signature of optically thick clumps, and t^/t^ ~ 1.5 h as been shown to be co nsistent 
with calculations using the analytic vorosity models by 'Sundq vist et al.l ( 2011 b'). 



The fact that clumps seem to be optically thick in UV resonance lines may then 
also help explain the unexpected weakness of observed UV lines such as Pv (Sect. 14711 ). 
as demonstrat ed by first a ttempts to in clude corresponding ef fects in line di agnos- 



tics dOskinova et al.ll2007l : iHillied l2008l : iSundqvist et al.. .201Q. 1201 Ibl : ISurlanI 12011) . 



Sundqvist et all (1201 Ibl) analyzed Pv and Hq. in the Galactic 06 supergiant A Cep, ac- 
counting for structures in density as well as for a non-monotonic velocity field. As is 
rather unaffected by optically thick clumpi ng in this O star, where as the strengths of the 



Pv lines are substantially reduced (see also lOskinova et aLll2007L for similar results for 



( Pup), the derived mass-loss rate is now only a factor of two lower than the prediction 
bv lVink et al.1 (120001) . This significantly alleviates previous discrepancies, and suggests 
that only mild reductions of current theoretical mass-loss rates are necessary. 

However, while the middle panel of Fig. |3] demonstrates that vorosity is indeed 
present also in LDI simulations, the resulting absorptio n redu ction is generally less 
than needed to fully explain the Pv observations (see also lOwoc ki 2008). This may in- 
dicate that the present generation of LDI simulations do not properly resolve the internal 
clump velocity structures, and thus overpredict the clump velocity spans 6v. Moreover, 
another problem with these simulations is that they do not show enough struc ture in 



the inner wind to reproduce the observed Hq, emission (ISundqvist et all 1201 Ibl) . We 
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conclude this review by briefly discussing this discrepancy between predicted structure 
and observationally inferred clumping factors in the inner wind. 




Figure 3. Left, middle: X-ray (left) and UV resonance (middle) line profiles for 
optical depths typical for ^ Pup at ~ 15 A and for the Pv doublet. Adapted from 
ISundqvist et al.l (l2011albl). R ight: Radial stratification of the clumping factor as 
derived for ^ Pup bv lNaiarro e t al. ( 2011 ) (black) vs. theoretical predictions averaged 
over ~ 1000 1-D LDl snapshots (blue). The same 1-D LDl simulations and synthesis 
calculation methods as in Fig.[T]are used. 



5. Clumping in the inner wind 

The right panel of Fig. [3] compares predictions from LDl simulations with the ra- 
dial stratification of fr] for ^ Pup , as rec ently derived from the comprehensive multi- 



diagnostic study by iNajarro et al.l ( 201lh . While the agreement is good both in terms 



of general behavior and in terms of absolute amount of clumping in the outer wind, the 
observationally inferred clumping law peaks at much lower radii than the theoretical 
one. 

Considering the multitude of independent observational studies indicating similar 
results (Eversberg et al. 1998; Bouret et al. 2003; Puis et al. 2006; Cohen et al. 2011), 
this motivates a re-investigation of the onset of wind structure predicted by standard 
LDl simulations. Note that while the LDl simulation displayed here does include base 
perturbations, it still relies on the SSF formalism for calculating the line force and 
the associated diffuse line drag (Sect. |2]). This probably overestimates the damping of 
perturbations close to the wind base, as suggested by first attempts to incorporate a more 



sophi sticated escape integral source function (EISF) into LDl models (O wocki & Puis 



Il999l) . Along with suitable base perturbations, future work will implement various 
forms of such EISFs into a new generation of LDl simulations, and examine to what 
extent this may induce substantial structure also near the wind base. 

Acknowledgments. J.O.S. gratefully acknowledge funding from NASA ATP grant 
NNX11AC40G 
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